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Abstract 

We discuss some of the difficulties that have been mentioned in the liter- 
ature in connection with the Bethe ansatz for the six-vertex model and XXZ 
chain, and for the eight-vertex model. In particular we discuss the "beyond 
the equator" , infinite momenta and exact complete string problems. We show 
how they can be overcome and conclude that the coordinate Bethe ansatz 
does indeed give a complete set of states, as expected. 
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1 Introduction 

There are proofs in the literature of the combinatorial completeness of the 
Bethe ansatz|l| for the six- vertex model and XXZ chain: |2|, ^ i.e. that for a 
lattice of N columns it gives all 2^ eigenvectors (states) of the transfer matrix. 
In the presence of an arbitrary field H there seems to be no doubt that this 
is so, and this is in agreement with the results of the numerical experiments 
we report in section 2. Other studies have been made, also indicating the 
completeness of the Bethe ansatz. [Q, ^ 0, |8| 

Even so, there still appear papers that either question this completeness, 
or at least appear to question it, when the field is zero, or at special values of 
the crossing parameter A (or r]). Statements have been made that "the Bethe 
vector vanishes" for states with more down arrows than up arrows, i, in and 
that it is incomplete or "singular" if some of the momenta are infinite. |0, 



12, 13| Here we show that these problems can be overcome in the coordinate 
Bethe ansatz if one recognizes that one is dealing with a set of algebraic 
equations, so must include an appropriately generalized "point at infinity" in 
one's considerations, and properly normalize the eigenvector. 

Recently it has been claimed that "Bethe's equation is incomplete" at 



special "roots of unity" values of A. [14] - [0] By this it is meant that the 
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Bethe zeros vi, . . . ,Vn, or equivalently the Bethe momenta fci, . . . , are not 
uniqely defined. They contain at least one exact complete string, and one is 
free to choose each string centre at will.Q This freedom is noted in |17], after 
equation (1.35) therein. We show that this freedom is because the eigenvalue is 
degenerate, so the eigenvector itself is not uniqely defined. Any allowed choice 
of vi,. . . ,Vn gives a valid eigenvector. The set of such choices is a curve in 
the eigenspace. For the simplest case, which is when vi, . . . ,Vn form just one 
single string, we show that the vectors on this curve span the eigenspace. Thus 
the Bethe ansatz is complete for this case, precisely because of this lack of 
uniqueness. We fully expect this argument to generalize to more complicated 
cases. 

In all the cases we have looked at, we have found that the Bethe ansatz is 
indeed complete: it gives all the eigenvectors (states). More precisely, it can 
be used to construct a basis for each eigenspace. 

Apart from the difficulty mentioned in section 6, the problems we en- 
counter can be resolved by the methods mentioned after equation ( ^6| ) and 
reviewed in the summary. 

We also present the coordinate Bethe ansatz equations for the eight vertex 
model in zero field, with an even number of columns, and discuss how the 
infinite momenta and exact complete string problems can be resolved. We 
expect these equations to be similarly complete. 

We further show that the functional relation between the eigenvalues 
T(v),Q{v) of the T and Q matrices can itself be written as a generalized 
eigenvalue problem, with the Fourier coefficients of Q{v) being the elements 
of the eigenvector. 

In one sense these problems have to be resolvable, since if the ansatz is 
complete for arbitrary field H and crossing parameter A, then (at least in 
principle) one can always deal with difficult cases by taking a limit. (When A 
is real or pure imaginary, and the field H is zero or pure imaginary, then one 
can choose the spectral variable v so that the transfer matrix is hermitian, so 
we know it is diagonalizable. Since the eigenvectors are independent of v, this 
means that it is diagonalizable for all u.) However, the question is whether one 
can do better than this and find all the eigenvectors for any particular values 
of H and A. For the problems that have been addressed in the literature, we 
believe the answer to be yes. 

To paraphrase Mark Kac, this paper is intended to fill a much-needed gap 
in the literature. 



-'^ Except that if = 0, iV is even and n = N/2, then one can require that the eigenvector g also 
be an eigenvector of the arrow reversal operator R. If one does this, then wi, . . . , w„ must satisfy 
the constraint (151). 

^This contradicts the statement in the abstract of |jl^, and repeated in the introduction of [Q, 
that "the Bethe ansatz equations determine only the eigenvectors which are the highest weights of 
the infinite dimensional SI2 loop algebra" . 
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2 The six-vertex model in a field 



The six-vertex model |19| - [27| has Boltzmann weights uJi, . . . ,ujq. Let 





Let cr = fJi, ... , cJiv denote the state of a row of N vertical arrows (+1 for 
an up arrow, —1 for a down arrow). Then the transfer matrix for a lattice of 
rows is the 2^ by 2^ matrix T with elements 

T^^a' = Trace i?^^ . (2) 

Considered as a function of wi, . . . , wg, it has the symmetries: 

T(tJi, . . . ,cj6) = (-l)^r(-wi, . . . , -We) = r*(u;4,u;3, u;2,t^i,t^6, ws) , (3) 

the superfix t denoting transposition. 

In statistical mechanics one wants to calculate the partition function 

Z = Trace T^"- , (4) 

whrer Mr is the number of rows of the lattice. It is therefore desirable to 
diagonalize the matrix T. This problem has been solved (for general values 
of u>i, . . . ,a;e ) by the Bethe ansatz |23| - [p6|.Pl 



In the Bethe ansatz one characterizes the state cJi, . . . , (Jat by the positions 
X = xi, . . . ,Xn of the down arrows, i.e. aj = —1 iff one of xi, . . . ,Xn is equal 
to j, else aj = +1. Because of the "ice rule" that there be two arrows into 
each vertex, and two arrows out, the number n is conserved, being the same 
on all rows of the lattice. For a lattice of N columns, the transfer matrix T 
therefore breaks up into + 1 diagonal blocks, one for each value of n from 
to N. Within block Ti, an eigenvector g has elements ^(^) — q{^x\^ • • • ? -^n) 
for each state a or X. 

The Bethe ansatz is the following guess at the elements of the eigenvector: 

g{xi,. ..,Xn) = Y. ^(^'i' ■■■^Pn) e'''^^''' ■ ■ ■ e'^^-''" (5) 
p 

where the sum is over all the n! permutations {pi, . . - Pn} of {1, . . . and 
the integers xi, . . . , x„ lie in the range 

1 < Xl < X2 < • • • < x„ < iV . (6) 

Substituting this ansatz directly into the eigenvalue/eigenvector equations, 
one finds the following sufficient conditions for g to be an eigenvector: 

) = , (7) 



■^In fact, ||2q| is more general yet, considering an inhomogeneous model where the field H and 
the rapidity variable v vary from column to column. Here we shall only consider the homogeneous 
model. 
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JNk 



PI 



A{p2,. . . ,Pn,Pl) = A{pi,. . . ,Pn) , (8) 



for j = 1, . . . , n — 1 and all permutations {pi, . . . ,Pn}- 
Here 



Si. 



We emphasize that (||) - @ are sufficient conditions (together with g 0) 
for g to be an eigenvector. They involve the Boltzmann weights only via the 
ratios 

(wiW2 -a;5C<J6)/(i^iW3) and c<j2^4/(^i^3) • (10) 

This implies that the transfer matrices of two models with different weights, 
but the same values of these ratios, commute. This can be proved directly by 



appropriately extending the method of section 9.6 of |2?], or of |^ or [29|). 
The eigenvalue corresponding to this eigenvector is 



A — LV-^ II 7 jr-\ r ^^4 I I / ik \ 

(11) 

If the weights wi, . . . , are all non-zero, it is convenient to write them as 

H+V -H-V H-Vi, 

uji = e a , u;2 = e a , UJ3 = e 

a;4 = e~^~^^b , uj^ = uoq = c . (12) 

Here H \s a, dimensionless electric field in the horizontal direction; V is the 
field in the vertical direction 
Defining 

A = (a2 + 62 _ c2)/(2a6) , e'^^ = e-^^ e'^^ , (13) 

and dividing Sjm by a constant factor e^^ab that cancels out of (0), the above 
equations involving cji . . . , simplify to 

Sj^rn = 1 - 2Ae*'"'" + e*(''^+'"'") , (14) 



A = e^N-2n)v I ^NH TT ah + {c^ - b^)e^-^ A a'-^-ahe'-^ 

(15) 

Note that V enters the above equations only via a factor 2")^ the 
eigenvalue A. This is because of the ice rule: T commutes with the diagonal 
matrix Sz which has entries (cJi + (T2 + ■ ■ ■ + o"Ar)5o-,o-'- The effect of V in the 
Boltzmann weights (^ is to pre- and post-multiply T by exp(S'2y/2). Within 
the diagonal block n, the diagonal entries of Sz are — 2n, so the effect of 



■*The weights ^5,0^6 are sinks and sources of horizontal arrows, so only occur in the partition 
function and transfer matrix in the product combination lo^^loq. This means that there is no loss 
of generality in choosing ujf, — oj^, — c. 
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introducing V is to merely multiply all entries of T by g^^"^")^. Without 
loss of generality we shall from now on take V = 0. 

The next step is to derive what Fabricius and McCoy call "Bethe's equa- 



tion". |15, eqn. 1.2] This is easily done when the sj^m sue all non-zero, but 
much of the apparent confusion in the literature arises when some of them are 
zero, so we proceed carefully. To be an eigenvector, g must be non-zero, so at 
least one of the coefficients A{pi, . . . ,pn) must be non-zero. In Appendix A 
we show that this implies that 



e'^'^ n ^i.™ = (-1)""' n ^--i (16) 

for j = 1, . . . , n. Hence these n equations are a necessary consequence of the 
linear equations (0) and ^ for the n! coefficients A{pi, . . . ,Pn)- We discuss 
below their sufficiency for the case when all of the sj^m are non-zero. 

Also, substituting all n cyclic permutations of {pi,... into (|8|), one 
gets n linear homogeneous equations for the n coefficients A that occur. Since 
at least one of them is non-zero (for some permutation {pi, . . . ,pn}), the 
determinant of coefficients must vanish, giving 

^iN(ki+k2-\ \-k„) _ ^ ^-j^y-j 

We refer to (|5|) - (|l7|) as "the Bethe ansatz equations" . 



Transformation to difference variables 

Define p, A, v so that 



a = psinh[(A — u)/2] 

b = psinh[(A + w)/2] (18) 
c = p sinh A 



then 



A = - cosh A (19) 
Define also vi, . . . ,Vn so that 

gifc,- ^ g2Hgi«,. ^ ^2H - _ /20) 

Then 



sinh A sinh[(z;j — vj + 2A)/2] 
" sinh[(wi + A)/2] smh[{vj + A)/2] ' 

If we define functions (j){v), Q{v) by 

^{v) = p^smli^{v/2) , 



(21) 



Q{v) = l{smh[iv-vj)/2] , (22) 
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then (15) can be written in the form 



_ , , e^^ 0(A - v) Q{v + 2A) + e"^^ ^{X + v) Q{v - 2A) 

^ - ^ Q(^;^ • ^^^^ 

We shall also introduce a parameter g by 

9 = -e^ , A = (<? + g-i)/2 . (24) 



This g is the -q of [15|. 



Commuting transfer matrices 

In the notation (jl^) - (|^), the ratios (^) that enter the eigenvector calcula- 
tion are 

2e-2^A and e"^^ . 

Keeping p, A, H fixed, and writing the transfer matrix T as a function T{v) 
oi V, it follows that 

T{v)T{v') = T{v')T{v) (25) 

for all v^v' . The transfer matrices T{v), T{v') commute. 
Further, if we define the Pauli matrices 

acting on the spin (arrow) in position j, then the logaritmic derivative of T[v) 
ai V = —A is a linear combination of the identity operator and the hamiltonian 

1 ^ 

^ = -o E{^oshF (cTjcTj+i+cT|cT|^i)-i sinh/7 (ctJctJ^, - ajcxj+i) + A } 
i=i 

interpreting suffixes + 1 as 1. Hence T{v) also commutes with this hamil- 
tonian. 

Incrementing v by 27ri is the same as negating all of a, 6, c.0 This merely 
multiplies T by (—1)^. The lowest and highest powers of that can occur 
in the expansion of an element of T{v) are —N and N . It follows that T{v) 
can be expanded in the form 

N 

T{v) = Y^Tre^^-^-)-"/^ , (27) 

r=0 

the coefficients T^- being matrices independent of v. 

The commutation relations (|25| ) imply that Tq , . . . , Tjq all commute with 
one another, and with Ti.. 



^Remeniber that and ujq always occur in pairs, so negating c has no effect on T . 
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Hermiticity of T and unitarity of P 

Negating v and H interchanges uji with 014, and 102 with ^3. From this 
transposes the transfer matrix T. If A, p are real and H, v are pure imaginary, 
or if V is real and p, H, A are pure imaginary, this implies that T is hermitian. 
Thus it is diagonalizable. Combining this with the commuation properties 
above, it follows that there exists a unitary eigenvector matrix P such that 

P^TrP= and P^HP = diagonal . (28) 

In both cases A is real. For these cases, the eigenvalue A must be real. 
Since H,v/X are pure imaginary in both, (p3|) suggests that Q{v + 2A)* = 



Q{v — 2A). From (22), this implies that: 

(i) A real: vi, . . . ,Vn are either pure imaginary or occur in pairs Vj, —v* 
positioned symmetrically about the imaginary axis; 

(i) A pure imaginary: vi,...,Vn are either real or occur in complex 
conjugate pair. 

In both cases it follows that the wave numbers ki, . . . ,kn are either real 
or occur in complex conjugate pairs. Similarly for ki, . . . , «;„. 

Note that P is independent of v, so provided H is pure imaginary, we have 
provedthat T = T{v) is diagonalizable and P is unitary for all complex v. 

Of course, we usually take the Boltzmann weights wi, . . . , to be positive 
real, so from that point of view we would like to take H to be real, rather 
than pure imaginary. However, it is reassuring to have a proof that T is 
diagonalizable even if the proof only holds when H is pure imaginary. This 
does include the "zero- field" case of main interest, as well as other problems 
thta have been looked at, such as the critical Potts model with q < 4.^ 

This implies that T will be diagonalizable for all H except posssibly for 
isolated non-zero values off the imaginary axis. Our question is: does the 
Bethe ansatz give all the eigenvectors? 

The matrix Q{v) 

Each eigenvalue A must have an expansion corresponding to (|^) : 

N 

A{v) = tre^''-^'>/^ (29) 

r=0 

Thus A is an entire function of v and the RHS of (^) must vanish when the 
denominator does, giving 

e^^0(A - Vj) Q{vj + 2A) + e-^^0(A + vj) Q{vj - 2A) = , (30) 

for j = 1, . . . , n. 



These are precisely the equations (16). The author noted (for H=0) in 



1971 that they imply the existence of a matrix Q{v) with eigenvalues Q{v) 
that commutes with T{v) and satisfies the matrix functional relation 

T{v)Q{v) = e^^(P{X-v)Q{v + 2X') + e-^"(l)iX + v)Q{v-2X') , (31) 



^In section 6.2 of |Q it is shown that this model is equivalent to a six-vertex model with 
a boundary seam. This seam is equivalent to introducing a horizontal field H — 9/N, where 
= 2 cosh0, so 9 and H are pure imaginary when q < 4 and the model is critical. 
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where A' = A — itt. This proved to be the key to solving the eight-vertex 
model. PI, M 



For the zero-field eight- vertex model an explicit construction for Q{v) is 
given in section 6 of [^] and in section (10.5) of |2?].[| This is specialized to 



the six- vertex model in equations (8), (96) and (97) of [Sjl, for H = and N 
even. 

The vectors $„(f |(t) therein are generalizations of the special eigenvectors 
^p we discuss in section 3. They form the columns of a matrix $„(t;). Writing 
[iV, m] = N] / {m\{N — m)\) , this matrix has [A'^, n] rows and [N,N/2] columns. 
Provided the rows are linearly independent, we can define the [N, n] by [N, n] 
matrix Q{v) by 

Q{v)^nivo\a) = ^niv\a) (32) 

vq being an arbirary fixed parameter, for all n from to N. Then it is shown 
in 0) and H that 

Q{v)T{v) = T{v)Q{v) , (33) 

so we can simultaneously diagonalize both T{v) and Q{v). Doing this, we find 
(for H = and even) that the eigenvalues Q{v) must indeed have the form 
given in equation (|4^) below. 

In the cases discussed in section 5, where some of vi,...,Vn form one 
or more complete strings, the eigenvalue A of T(v) is degenerate. This is 
reflected in the fact that the string centres (the average value of the Vj within 
a string) are not determined by the Bethe ansatz. One of the main objects of 
this paper is to stress that this is not a deficiency or "incompleteness" of the 
Bethe ansatz, but rather a strength. It means that it can be used to construct 
a complete basis of the eigenspace. 

However, the corresponding eigenvalues of Q{v) are not degenerate, so one 
can fix the string centres by requiring that g in (^) be also an eigenvector of 
Q{v). This is what Fabricius and McCoy have achieved. In physicist's terms, 
they have resolved the degeneracy of the eigenvalue A; in mathematician's 
terms, they have made a particular choice of the basis of the eigenspace of 

nv). 

Beyond the equator: the relation between the n and N — n 
solutions 

As we note below, there appears to be no problem solving these equations 
in the presence of a non-zero field H (pure imaginary or real), even when 
n > N/2 and there are more down arrows than up. However, the eignvalues 
A are then the same as the mirror case (where all arrows are reversed), with 
n ^ N — n and H — > —H. Is there a relation between the two solutions? 

This problem has been studied by Bazhanov, Lukyanov and Zamolod- 
chikov p^, p3|. Consider the equation (23) for some n and H, with a function 



Qi{v) = Q{v), and the same equation with n, H replaced hy N — n, —H, with 
a different function Q2{v) = Q{v) but the same A. Eliminate A between the 



we remark in section 7, there is a typing error in eqn (10.5.8) of [g^- '^j+i therein should 
be CTj-i. 
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two equations. We obtain 

W{v + \) = W{v-\) , (34) 

where the "Wronskian" W{v) is defined by 

. ^ e^^Qi(z; + X)Q2{v - A) - {-l)^e-^^Q,{v - X)Q2{v + A) 

(l){v) 

(35) 

We continue to require that Qi{v) and Q2{v) be of the form (^), with n 
replaced by n, — n, respectively. It follows that 

W{v + 2TTi) = W{v) . (36) 

For arbitrary A, excluding the "root of unity" cases discussed in section 5 
in which iX is a rational fraction of vr, the only solution of both ( |3^ ) and ( |36|) 
is that VF(f) be a constant D. Hence for non-"root of unity" cases, Qi{v) 
(with parameters n, H), and Q2iv) (with parameters N — n, —H) must satisfy 
the relation 

e^''Qiiv + X)Q2{v-X)-{-l)^e-^''Qiiv-X)Q2{v + X) = D,p{v) . (37) 

For H = 0, this is equation (47) of (0]). We shall discuss it further below, 
particularly for the case when H = and is even. 



Pronko and Stroganov|10] have also addressed this problem, in a different 
manner. Using their terminology, we shall show in section 4 that it is indeed 
interesting, and perfectly possible, to consider Bethe's equations 'on the wrong 
side of the equator'. 

Bethe's equation as a generalized eigenvalue problem 

We can write (|23) as 



A{v)Q{v) = (-l)"{e^^(A(A - v)Q{v + 2A) + e-^^(^(A + v)Q{v - 2A)} (38) 
and (p{v), Q{y) as 

N 

m = Y: /.e(^-2'-)-/2 (39) 

r=0 

n 

Q{v) = Y,q,e^^-^^>l^ (40) 
j=0 



Substituting these expansions, together with (|29| ), into the A,Q equation 
above, one obtains 

n 

E^m9j=0' (41) 

3=0 

where if < i — j < A^, 
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(42) 

else bij = 0. Here i = 0, . . . , N + n. 

We see that we have N + n + 1 equations for the N + n + 1 unknowns 
to,...tN and qo : qi : • • • : Qn- These are an alternative form of Bethe's 
equations. They define the eigenvalue A{v) and (usually) the function Q{v). 
They are linear in the tj, and homogeneous and linear in the qj. The tj play 
the role of a set of "eigenvalues" , qj that of an "eigenvector" . 

This form of Bethe's equation has some advantages which we shall mention 
as we come to them in the following four sections. In particular, suppose that 
we actually know, or have guessed, the eigenvalue A(f ), and hence to, ... ,1^. 
Then the equations are a set of homogeneous linear equations for qo, . . . ,qN, 
and can be solved by the standard apparatus of linear algebra. Let B be the 
N + n + 1 byn + 1 matrix with elements bij . Then we can distinguish three 
cases: 

1) B has rank n + 1: then there are no solutions for qo, . . . , qn- A(f ) is not 
an eigenvalue. 

2) B has rank n: there is one solution for the ratios qo : qi : ■ ■ ■ : qn- This 
presumably means that the eigenvector g is unique: then A{v) is an eigenvalue 
with degeneracy one. 

3) B has rank less than n: there is more than one solution for qQ : qi : ■ ■ ■ : 
qn- The eigenvector g is presumably not unique: then A{v) is an eigenvalue 
with degeneracy greater than one. 

Thus we can use these simple considerations to determine whether a given 
eigenvalue is single or multiple. 

"Bethe's equations" 

There are a huge number of equations in (|^): (n — 1) x nl/2 homogeneous 
linear equations for the n! coefficients A{pi, . . . ,Pn). Fortunately it seems that 
they always permit at least one (possible more - this is the source of some of 
the misunderstandings in the literature) non-identically zero solution. This is 
of the form 

A{pi,...,pn) = epC-' n ' (43) 

l<i<j<n 

where ep = ±1 is the sign of the permutation and the tij must satisfy 

tij Sji — tji Sij . (^^) 

At least one of tij and tji must be non-zero, else the A(j)i, . . . ,pn) would all 
vanish and g would be the zero vector. 

If the Sij are all finite and non-zero and the kj are finite, then we can take 
tij = Sij and choose the normalization factor C to be unity. The problems 
discussed in this paper arise when this is not so. (Apart from the equal Vj 
difficulty touched on in section 6.) For these cases one should choose C so 
that the maximum term in the summand in (|5|) is finite and non-zero (say 
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unity). This maximum is to be taken over all permutations P and all values 
of Xi, . . . Xn allowed by (n). 



This is the solution given in (8.4.10) of |27|, except that there we took 



tij = Sij for all We wish to be more general here so as to cope with the 
situation when some of the sa vanish. 



Substituting (|4^) into we obtain the n equations 



( 1) im,j /ij,r 

m=l,mj^j 



(45) 



for j = 1, . . . ,n. Together with (p4|), this implies (16). 
We remind the reader of equation (14), namely 



1 — 2Ae^^'~'^^ + Q^'y^i+^j)"'^^ 



(46) 



We refer to (E3) - (E^) as "Bethe's equations" , which is a slight extension 



of the terminology of Fabricius and McCoy |15, eqn. 1.2]. They are sufficient 
conditions for (|^) and (|8|) to have a non-zero solution for the coefficients 
A{jpi, . . . ,Pn)- They form a set of coupled equations for the e^^ 



ratios tij '. tji- 



'J . Sij and the 



Apart from section 6, the problems we shall be discussing occur when some 
of these variables are zero or infinite. The resolution is always to re-express (|5|) 
and (43) - (|45| ) in terms of finite combinations of powers of the variables e^^^ 



and tij, to solve the equations simultaneously for these, and to allow for the 
possibility of a solution containing one or more arbitrary degrees of freedom. 



2.0.1 Momentum 

When V = —A, e^'^"'~^^^T{v)/c^ is the matrix with entries Hi '^('^ii <^i+i)- It 
shifts all arrows one column to the right. Doing this to the eigenvector g is 
equivalent to multiplying g by exp[i{ki + • • • + kn)]. Similarly, the matrix 
e(^-2'^)-^r(A)/c^ shifts all arrows one to the left. Hence, writing A as A(v), 

and e*^'^^"' ^^"'^ must be an A^th root of unity, in agreement with ([T7|) . From 
(p^) it follows that 

gi(fci+-+fcn) ^ e2"^Q(A)/Q(-A) . (47) 



Numerical calculations 

Method 

To fix our ideas, we conducted a number of numerical experiments on the 
above equations for lattices of small size: up to A = 8. We fixed H, A 
and p and evaluated the matrix coefficients Tr in (p7|). We then assigned v an 
arbitrary value and diagonalized T{v) directly to obtain the eigenvector matrix 
P, and verified that it did indeed diagonalize all the T^. For each eigenvalue 
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A this gave us the coefficients t^. in (|29|). We constructed the matrix B and 
determined its null space, and hence all solutions of (^T|) for the qj. This 
gave us the function Q{v). We calculated its zeros to obtain the vj from (p2|). 
We then calculated the kj and the wave numbers kj from (pO[), and then the 
Sij from (46). Usually the Sij were all non-zero and we took tij = Sij, and 
calculated the coefficients A{pi, . . . ,pn) from ( |4^ ) and the elements of g from 
(^). Finally we normalized this vector (so its largest element was one) and 
compared it with the correspondingly normalized column of the matrix P. 



Results for i7 ^ 

We first took H to be non-zero, either real or pure imaginary, and A to also be 
either real or pure imaginary. At this stage we avoided the "roots of unity" 
cases when iX is a rational fraction of vr. We encountered no problems with 
the above procedure. For a given number n of down arrows, we found no 
eigenvalues of T{v) that were identically degenerate for all v.^ The column 
nullity of B was always 1, so there was only one solution (to within normal- 
ization) of (|4^) for the qj. The Sij were all non-zero and the wave numbers ki 
were all distinct so that ^ gave a unique non-zero eigenvector. We worked to 
about 17 decimal digits of accuracy, and the error in the eigenvector elements 
was no bigger than 10~^^. 

For H pure imaginary, we also observed that the eigenvector matrix P was 
unitary, that for every Vj there was a conjugate according to the rule given 
after eqn. ([2^), and that the wave numbers ki,...,kn were either real or 
occurred in complex conjugate pairs. The functions Qi{v), Q2{v) had n,N — n 
finite zeros respectively, and satisfied (^). 

3 Particular values of H^X: some very spe- 
cial eigenvectors 

Here we consider the case when H, A satisfy the relation 

(_eA±2/^^^ = 1 , (48) 

and show that the Bethe ansatz then admits some some very special eigen- 
vectors. They are the analogues of the special eigenvectors of the zero-field 
eight-vertex model obtained by the author in section 7 of [pl|]. 

Almost the simplest ansatz one can imagine for an eigenvector ■0 of the 
2^ by 2^ transfer matrix T{v) is the direct product form 




where gi, . . . ,g]\f are some parameters to be determined. 

^Obviously they can be degenerate for special values of v. if v — ±A, then c^^T{v) is 
"momentum" operator that shifts all arrows in a row one column to the right (or one to the 
and has eigenvalues which are iVth roots of unity. 
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Define q by (|24[): q = —e, and let 

q = q^^ , (50) 

making one of the two possible sign choices here and in the following equations. 

Following the method of section 9.8 of 1 27], we find that T{v)ip has a simple 
structure if 

9j = {e^''qyg for j = l,...,iV , (51) 
except that we need the cyclic boundary condition gN+i = 9i' this implies 

{e'^^qf = 1 , (52) 
which is (48). The parameter g is arbitrary: we can choose it at will. Hence 



we can regard ■0 as a function ■0(5) of g. 
Then we find that 

Tiv)i^{g) = u;^i;{g')+u;Uig") , (53) 

where 

g' = Q9 , g" = Q'^g ■ (54) 

Now look at the sub-space with n down arrows and Sz = N — 2n. Then 
i^ig) — g^4'n, where, in terms of the positions xi, . . . ,Xn of the down arrows, 
0n is a vector with entries 

</>„(xi,...,x„) = (g)(-i+-+-") . (55) 

In this sub-space ( |53| ) becomes 

r(T;)</.„ = (a^e^^g" + 6^e-^^r")<An • (56) 

Thus 6 n is an eigenvector of T(^v), with eigenvalue 

A = a^e^^g" + 6^e-^^r" • (57) 



This is true for all n from to N, provided only that the restriction (|48|), or 
more specifically (|5^, is satisfied. 



Reconciliation with the Bethe ansatz 

How can we reconcile this with the Bethe ansatz? Simply by taking 

Vj ^ ^oo for j = 1, . . . , re . (58) 



Then, from (gg), 

gifc, ^ g2H~ fQj, j = l,...,n , (59) 

so all the exponential factors containing xi, . . . ,x„ in (^) are proportional to 
(pnixi, ...,Xn) and 

g (X(j)n ■ (60) 



Also, from (|2^) and (^), we obtain the result (|57D for the eigenvalue A. 
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It would seem that there is nothing more to say: the vector is a special 
case of the Bethe ansatz when all the Vj tend to =Foo. The eigenvalue is given 
by (HD. 

However, the alert reader will notice that we have said nothing about 
Bethe's equations. More seriously, all the kj are finite and equal. This is a 
potential problem in itself, since a casual inspection of (|4^ suggests that it 
implies that tj j = tj^i for all i,j from 1 to n. If we then take A[P) to be 
given by (^) and substitute into (^), all terms will be equal except for the 
sign factor ep. They will therefore all cancel and we shall obtain g = Q. (We 
return to the general problem of what happens when fcj = fej, Sjj- = Sj^i / 
for some pair of values in section 6.) 

However, e**^^ = q, so from ([l^ and (|2^ , 

Sj,m = for all j,m . (61) 
At first sight this appears to only make matters worse. If we use the 



"normal" solution of (f44D, namely tj^m = sj^m, then (^) gives A{P) = 0. 
Now each term in (^) vanishes, not just their sum! 

The answer is of course that one should not use this solution. In fact there 
is now no reason to use the subsidiary "ansatz" ( |43| ) at all. For any finite 
choice of coefficients A{P), the equations (0) are satisfied simply because the 
Sj,rn vanish. All that remains is to satisfy (|^. One simple choice that does 
this is to take 

A{P) = 1 (62) 



for all permutations P. From (|5^ ) and (pg|), 

e'^''^ = 1 , (63) 

so both (^) and (^ are satisfied and (^) becomes 

g{xu...Xn) = n!e^'=(^i+-+^") , (64) 

where k = ki = ■ ■ ■ = kn- Apart from the non-zero normalization factor n! , 
this is the result (Issj) . 

These special eigenvectors are a good illustration of how one can satisfy 
the original Bethe ansatz equations when some of the Sjj vanish. The form 



(42) is not part of the original ansatz, but an addition to it. It is a necessary 
consequence of if all of the Sij are non-zero, but if enough of them vanish 
it ceases to be necessary. 

Still, if ( |43| ) is true in general we should not be too ready to abandon it 
in particular. We can still take A(P) to be given by (^3|). Equation (|4^ ) now 
imposes no restriction on the tij, but (^) does. For the case that we are 
discussing in this section, a simple solution of ( |45| ) that avoids the problems 
that occur when tij = tj^i is to take tij = —tj^i ^ for all i,j. Then A{P) is 
independent of P and we regain the solution just discussed. 

Another solution can be obtained, under the more specialized conditions 
(p^) and (p3[), by letting the Vj in section 5, for the case M = n, tend to =F oo. 

More generally, we can take the tij to be arbitrary and non-zero, with 
tij 7^ tjj,for 1 < i,j < n — 1, and then use (45) to determine the ratios 

ij,n/in,j • 
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In all three approaches, note that we use (^), rather than (pif), to deter- 
mine some of the ratios tij/tji. This is a basic feature of the algebra of the 
next two sections. 

Infinite VjS: the general situation. 

Suppose just r of the vi,. . . ,Vn equal — oo, s equal +00, and the remaining 
n — r — s are finite and arbitrary. More precisely: 

Vj = finite , for j = 1, . . . ,n — r — s , 

Vj = —00 , e***^ = q , for j = n — r — s + l,...,n — s , (65) 
Vj = +00 , e**^^ = q^^ , for j = n — s + 1, . . . ,n . (66) 

Let us call these three cases type 1 to type 3, respectively. From ([l^ and 
(p^), if Vj is of type 2 and Vm is not of type 2, then 

Sjm = {q- e'''^)/q , Smj = -q{q - e'^^) (67) 

so Smj/s /tjm = -q^- Similarly, if Vj is of type 3 and v^ is not of 

type 3, then Smj / Sjm — tmj/tjm — q ^• 

If Vj and Vm are both of type 2, or both of type 3, then Sjm = Smj = 
and A{P) is not necessarily given by (|4^). Strictly, we should go back to the 
original Bethe ansatz equations (0) and (^. 

However, for similar reasons to those above, it seems that we can without 
loss of generality take A{P) to be given by (p3|), so long as we realise that 
(^) no longer defines tjm/tmj for Vj and Vm both of type 2, or both of type 3. 
Then (I) gives the equation (|45|). Taking Vj to be of type 1, 2 or 3, we obtain 

^iNk, ^ g2s-2. j^n-r-s ^.Sm,j/s,,m) , j = I, ■ ■ .,n-r-S , 

e^^^q^= q'-~^^Ul=n-r-s+li-tm,j/hm) , J = U-T - S+1, . . . , U- S , 
e'^V^= q''-^''nm=n-s+li-tm,j/hm) , j = H- S + 1, . . . ,n , (68) 

where all three products exclude the value m = j. 

Taking the product of each of these three equations over the allowed values 
of j, the —Smj/sj^m, —tm,j/tj,m factors cancel out, leaving 

gi7V(feiH hA:„-r-s) _ ^2{s-r)(n-r-s) 

^2NrH ^ ^T{2n-2r-N) ^2NsH ^ ^s{N+2s-2n) ^gg^ 

Eliminating H gives 

^2rsiN+r+s-2n) ^ ^ ^ j^^^q^j 

SO we see that such infinite zeros can only occur when r,soTN + r + s — 2n 
vanishes, or when q (and e^^) is a root of unity. 

Let Q{v) be given by (|2^ , but with the product restricted to the finite Vj: 

n—r—s 

Q{v) = n sinh[(t; - v,)/2] , (71) 

3=1 
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then, to within factors independent of v, 

Q{v) = e^'-'>/^Q{v) . (72) 
Substituting into the eigenvalue equation (^), we obtain 

^ ^ ri)n ^ <t>{\ - v) Q{v + 2A) + u-^ 0(A + v) Q{v - 2A) 

Q{v) 



where 
Hence 



u; = e^"{-qY-' . (74) 



,2r _ r{2n-2s-N) 2s _ s{N+2r-2n) 



and we see that if r, s, + r + s — 2n are aU non-zero, then lo must be a root 
of unity. 

4 The zero-field model: H = 

The next step was to turn off the field, setting H = 0. For N odd no problems 
appeared: we were able to calculate all the eigenvectors in all the subspaces 
n = 0, . . . ,N without difficulty. For even we encountered two problems, 
both of which have been discussed previously in the literature. 

4.0.2 "Beyond the equator": even and n > N/2. Q{v) has 
infinite zeros 



In this "beyond the equator" case[|lO[, (|37D has the simple solution 

Q2iv)=Qiiv) . (75) 

Certainly (^Tj) permits this solution, and we observe numerically that B 
has nullity one, so it is the only solution. Hence Q{v) is the same for n as 
for N — n. This is quite consistent with (^) - it merely means that 2n — N 
of the zeros have gone off to infinity. This is even easier to see in (^0|): the 
first n — N/2 coefficients qj vanish, as do the last n — A^/2 coefficients. The 
degree of the Laurent polynomial is reduced from n to N — n, which is still 
consistent with (pl|) . 

From (pO|) and ([2^), this means that n—N/2 of the wave numbers ki, . . . ,kn 
are given by e*^-' = q, and another n — N/2 by e^^^ = 1/q. 

Faddeev and Takhtajan |p state that "Bethe's vector vanishes" for n > 
N/2. For n = 1 + A^/2 this is certainly not true of the vector g as given by 
(^) and (43). We have observed numerically that it is non-zero and that it is 



indeed the eigenvector of T. 

For n > 1 + A^/2 there is a problem, but it can be overcome, using the 
working at the end of the last section, taking r, s therein to be n — A^/2, so 
that (|69|) is satisfied for all X,q . When vj = Vm = ^oo, then Smj = sj^m = 0, 



so (44) tells us nothing about tm,j,tj^m- instead they should be chosen to 



satisfy the last two equations of (|68|). There will be many ways to do this. 
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corresponding to the fact that A{P) enters (P) only via its sum over ah ways 
of permuting VN-n+i, ■ ■ ■ , ■ ■ ■ > ^n- One simple way is to take 

tm,j = ~tj,m- The N — n finite zeros Vj are given by the first of the equations 
(|6^), which is the same as (45) when n is replaced by — n, i.e. the Bethe 
equations for the right side of the equator. The eigenvalue equations are the 
same for both n and N — n. 

The eigenvector equations are different, since we must include all n zeros 
in the product in (|4^). The resulting coefficients A{pi^ . . . ,pn) are finite and 
non-zero. The equations (^) and (^) are satisfied, so the vector g with elements 
(^) must be an eigenvector if it is non-zero. It is non-zero in our numerical 
experiments, and it is the eigenvector corresponding to the eigenvalue A. 

Thus despite the assertions that have been made in the past, the Bethe 
ansatz can be used to construct the eigenvector g for n > N/2, even when 
is even. It is a furphy that Bethe's ansatz does not work on the wrong side of 
the equator. 

4.0.3 even and 2 < n < N — 2: a single bound pair 

The other problem that we encountered first occurs for = 4 and n = 2, 
then for even A^ and 2 < n < N — 2. It is referred to by Bethe himself jl], after 



eqn. 23] and has been considered by others since |l3|],|3j, eqn.3.2.23b]. 
For some eigenvalues with momentum ±1, i.e ki + ■ ■ ■ + kn = or vr, we found 
that Q{v) had a pair of zeros vi,V2 such that vi = X, V2 = —A. From ( ^0|) 
this implies that 

More strongly, it was always true for such pairs that 

gi(fci+fe2) ^ _i ^ (77) 

so from (^6D, 

si2 = 2Ae-*'=i = oo (78) 

while S21 vanishes. 

Since si2, S21 do not both vanish, we can take as usual ti2 = S12, ^21 = S2i- 
In fact t2i vanishes strongly as this situation is approached (say by turning 

off the field H): from (||), 

«2i = C21 e^^^-^)'^! , (79) 



where C21 has a finite non-zero value given by (45). 

There is no problem solving Bethe's equations for the other k^, • • • ? 
and for C21. In principle one can then substitute these expressions into ( ^3| ) 
and (^) and extract the terms and elements that grow most rapidly as e*'^^ 
vanishes. Choosing C as stated after (|4^), each term in the summand of (|5|) 
will give a finite contribution to the eigenvector g. Some will be zero, but our 
numerical experiments indicate that the total vector g is not zero, and is in 
fact the correct eigenvector. 
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Numerical results 



For our numerical experiments, we simply assigned e^^^ the numerically small 
but non-zero value 10~^^i (the i being necessary for ki, k2 to be complex 
conjugates), calculated the Sij other than S21 from (|4^), then calculated S21 
from ( [45D , substituted the results into (^3|) and (|5|) and normalized the vector 
g. We obtained the correct eigenvector to approximately 14-digit accuracy. 

For = 4 there was just one such eigenvalue A, in the n = 2 central 
block. For = 6 there was one in the n = 2 block, two in the n = 3 block, 
and one in the n = 4 block. For = 8 there were 1,2,5,2,1 in the n = 
2, 3, 4, 5, 6 blocks, respectively. This suggests (tentatively) that the Catalan 
numbers may count such eigenvalues. The momenta were —1, except for a 
single eigenvalue with momentum +1 in each block with 3 < n < — 3. 

When N/2 + 1 < n < N — 1, there are eigenvalues where the two problems 
we have just discussed occur together, i.e. more than one of the e^^i are 
equal to g, more than one to and two of the others are itzcx). The above 
procedures were built into our computer program and worked perfectly, giving 
the correct non-zero eigenvector of T. 

For A^ even, we also calculated the matrix ^n{v) of equation (96) of |31]. 
We found that its rows were, as expected, linearly independent, so we were 
able to use ( |32[ ) above to calculate Q{v), and did indeed find that this matrix 
was diagonalized by the same matrix P that diagonalizes T{v), and that 
each eigenvalue was the function Q{v) discussed above (more precisely, the 
eigenvalue was Q{v)/Q{vq) ). 

The main lesson from this and the previous section is that Bethe's equa- 
tions should be viewed as a set of coupled non-linear equations for the e^^^ , s 



and the ratios tij : tji. We do not necessarily proceed by solving ( [46| ) and ([14[ ) 
for the ratios : tji. For some i and j it may be appropriate to obtain this 
ratio from (45). 



5 H = and q a root of unity: exact com- 
plete strings 

Now we come to the case discussed by Deguchi, Fabricius and McCoy in IQ 
- [0], where A is a rational fraction of iir, i.e. there exist integers u,M (with 
no common factors) such that 

A = iu7r/M , q = -e'^^/^ , (80) 

using dgl). 
Then 

= 1 , (81) 

and there is no smaller integer power of that equals one. 

For the moment we allow H to be arbitrary: we shall show that the further 
restriction ( |9l|) is necessary for a string (more precisely, a single string) to 
occur, and from then on take H to be zero. 
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The set of Bethe zeros vi, . . . ,Vn may now contain one or more "complete 
strings" , in which M of them, say vi, . . . ,vm, are related by 

Vj =vi + 2(j - 1)A for j = 1,...,M . (82) 

This implies that fj+i = Vj + 2A and vi = vm + 2A — 2ii^7r, so from ([2l|), 

Sl2 = S23 = • ■ ■ = SM-1,M = SM,l = . (83) 

We see that some of the Sij vanish, so we have to be careful with Bethe's 
equations. 

The set {vi, . . . , vm} may of course contain more than one complete string. 
For simplicity, from now on we shall restrict our attention to the case when 
there is only one string, but we fully expect our methods and comments to be 
applicable to the general case. 

For j = 1, . . . , M, the equation (^) becomes = 0, which is true, but not 
helpful. For j = M + 1, . . . ,n, both sides of the equation are non-zero. Using 
the form ( pl| ) of Sij, we obtain 

Sl,jS2,j ■ ■ ■ SM,j = (-l)^Si,lSi,2 • • • Sj^M , (84) 

SO for j > M, equation (p^) simplifies to 



e^Nk, ^ (_i)n-l-Af -Q s„,,,/s,,^ . (85) 
m=M+l,mj^j 

Taking the product of these n — M equations, the Sj^rn,Sm,j factors cancel, 
leaving 

^iN{kM+i+--+k„) _ _ |-ggj 

From (|l7| ) it follows that 

Now we note from ( ^2|) and (]2^ ) that, for j = 1, . . . , M, 

e^'^ = ^e^'^ ll'l:]"' , (88) 
1 + q'^^ 

where z\ = exp(i'i) and more generally 

= e}'^ = q^^-^zi . (89) 
Substituting this into ( |87[ ) and using (^), we obtain 

This is an extra condition on q and H that must be satisfied for a complete 
string to occur .0 From this and (81) we see that 



^At least for just one string to occur, but the same condition appears to be necessary for any 
number of strings. 
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An obvious and interesting solution of this equation is 



H = , (92) 

and from now on in this section we shall take H = 0, but we do note that 
there are other (pure imaginary) values of H for which complete strings may 
occur. 



If N is even, ( |90D is implied by (81). If is odd, the two equations 



together imply that q^^ = 1, but this is consistent with being a primitive 
Mth root of unity only if M is odd. Thus we have two possibilities: 

g^*^ = 1 , N even , 

q^ = I , iV and M both odd . (93) 
Apparent difficulties: (1) calculating g 

There are two problems that appear when a complete string occurs: the Bethe 
equations do not have a unique solution, and if we use the obvious solution 
tij = Sij of ( |4^ ) in (43), then every coefficient A{P) vanishes, so the eigenvector 



g, given by (^) also vanishes .|^ 

Let us dispose of the second difficulty first, since it is fairly straightforward. 
For the algebraic Bethe ansatz, it is considered by Fabricius and McCoy in 



the remarks after their equation (1.36) of |17]. 

We can still use the ansatz (|43| ) for the coefficients A{P) and attempt to 
satisfy the modified Bethe's equations (44), (|45|). Taking tij = Sij for all 
we note from (^3|) that 

^12 = *23 = • • • = tM,l = . 



For j = M + 1, . . . ,n, (|45|) becomes the above reduced equations 
which can be viewed as fixing vm+i, ■ ■ ■ ,Vn- 

Let vi be assigned arbitrarily. Then V2, ■ ■ ■ , vm and e*'^! , • • • , e*'^*^ are given 
by (ID and (H). 

We still have to satisfy (^5|) for j = 1, . . . , M. We can do this by using it to 
determine the ratios tj-ij/tjj+i, for j = 1, . . . , m (with to,i = tm,m+i = im,i)- 
We could of course have originally formulated Bethe's equations (which are 
just a set of algebraic equations and may well have solutions at zero or infinity) 
in terms of these non-zero, finite ratios. 

We also take C in (^) to be ti2. Then some of the coefficients A{P) will 
depend on ii2, ^23; • • • > tM,i only via their ratios, which are given by (|45|). The 
remaining coefficents A{P) will vanish. 

Thus all the coefficients are finite, some are non-zero, and we may hope 
that the eigenvector g, given by (|5|), will befinite and non-zero. In our numer- 
ical experiments this is what we have found. 

If there are 7 complete strings, then C should be the product of 7 factors 
tij, being one of the vanishing tij from each string. 



^"The only way the rhs of (^3|) could be non-zero would be for it to contain the factors 
S21, S32, ■ ■ ■ , sm,ai-i, si.M ■ But this cannot happen as there is no inverse permutation P' such 
that p'2 > p'itP'^ > P2, ■ ■ ■ ,p'i > p'm'- fhe inequalities are inconsistent. At least one of them must 
fail, which is the reason for the renormalization of A{P) proposed below. 
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Apparent difficulties: (2) non-uniqueness of the eigenvector 



We said above "let vi be assigned arbitrarily" . Why is this allowed? Should 
not its value be determined? This is the problem that concerned Deguchi, 
Fabricius and McCoy. - 

The answer is that for any value of vi the above procedure satisfies (^3|) 
- (p5|) and therefore the Bethe ansatz equations (^) - (p^. Provided g is not 
zero (and our numerical calculations indicate that it is not), then it must be 
an eigenvector of the transfer matrix T(v). We are free to choose vi as we 
wish. There is no a priory need to "complete" Bethe's equations.! 16 1 

Let us look more closely at what is happening. 

From (^) and (22), the M zeros vi, . . . ,vm contribute to Q{v) a factor 



M 



Yl sinh[(w - Vj)/2] cx sinh[M(t! - vi)/2\ 



(94) 



This factor cancels out of (p^), except only for a constant (— 1)*^. 

It follows at once that if Q(u) satisfies (p3|), and hence the Bethe equations, 
then so will any other function Q{v) with a different value of Vi. 

There is nothing remarkable about this. It means that the matrix B has 
rank less than n, so there is more than one solution qq, . . . ,qn of (|4l]). This 
in turn is a signal that the eigenvalue A{v) is degenerate, for all v. 

Of course, changing vi will change Q{v), so the matrix Q{v) will not 
be degenerate. If we construct it explicitly as above, then diagonalizing Q{v) 
rather than T(v) will resolve the degeneracies of T{v). Further, these will give 
the eigenvalues and eigenvectors obtained by taking the limit as A approaches 
the value (|80|). This is what Fabricius and McCoy have done, and the results 
are interesting. 

What we are concerned with here is showing that there is nothing in their 
work that indicates that the Bethe ansatz is incomplete, in the usual sense of 
not giving all the eigenvectors or "states". The fact that Q{v) is not uniquely 
defined by (^) - (17) is precisely the reason why one can use these equations 
to obtain a complete basis of the eigenspace of the eigenvalue A(u). 

A single string containing all the vi, . . . ,Vn- 



For simplicity we further restrict our attention to the case when all the 
vi, . . . ,Vn lie in the string, i.e. 

M = n . (95) 

However, we expect the substance of our remarks to generalize to n > M, and 
indeed to the case when vi, . . . ,Vn contain more than one string. 
The eigenvalue A is given immediately by (p^): 

A = g«[0(A-i;) + </>(A + ?;)] = q'^ia^ + h^) . (96) 

To within a sign, it is the eigenvalue for the "vacuum" state, when all spins are 
up. We shall now use the original Bethe ansatz equations (^) - ([T7|) to obtain 
explicit expressions for all the eigenvectors g corresponding to this eigenvalue. 
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for a given value of the number n of down arrows. We shall not use "Bethe's 
equations" (H3|) - ( |45|). 
From (^) and (^l|), 

Sl2 = S23 = ■ • • = Sn-l,n = Sn,l = . (97) 

The coefficients A{1, 2, . . . , n),A(2, 3, . . . , n, 1), . . . , ^(n, 1, 2, . . . , n — 1) en- 
ter the equations (0) only with multiplying factors Sij that belong to the set 
(p7|). This means that these coefficients do not enter at all. It appears that all 
other coefficients do enter with non-zero multiplying factors, and the equations 
ensure that they vanish. 

Choosing A{1,2, . . . ,n) = 1, from (|8|) it follows that 

^(i,j + l,...,n,l,...,j-l) = e^A'(fc.+fc,+i+-+M , (98) 
From (^) it follows that 

n 

g{xi, ...,Xn) = ^ exp{i[kiXn+2-j + k2Xn+3-j H 

J=l 

• • • + kj-iXn + kj{xi+N) + kj+i{x2 + N) + ■■■ + kn{Xn+l-j+N)]} . (99) 

The cyclic property ^(xi, . . . , Xn) = g{x2, ■ ■ ■ , Xn,xi + A^) is manifested by 

All of the Bethe ansatz equations (||) - (p^) are now satisfied. We still 
have the parameter vi, or equivalently zi, at our disposal. 

Substituting (^8|) into (|9^) and dividing by a common normalisation factor 
(1 — (— zi/g)"')/n(l + zi/q)^ , we obtain 

n n 

gizi\xi,...,xn) = n-i^g-^+-"+^(i-^')n(l + ^ir^^""'-"' , 

j=l r=l 

(100) 

taking Xn+i = xi + N and exhibiting the dependence of g on zi. 



From (93), g^" = 1, which means that the summand in ( |100| ) is unchanged 
by replacing j by j + n. Replacing j, zi by j — 1, q'^zi, we observe that 



g{x^zi\xi, . . . ,Xn) = g{zi\xi,. . . ,Xn) ■ (101) 



Because of the restrictions (|§), the RHS of (|100| ) is a polynomial in zi of 
degree N — n. Hence we can expand the vector g and its elements in powers 
of zi: 

N-n 

5(^1) = ^ Z^Ck , 
k=Q 

N-n 

g{zi\xi,...,Xn) = ^ z\ck{xi,...,Xn) , (102) 
fc=0 

Cfc being a vector with elements Cfc(xi, . . . , Xn)- 
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I have only verified this for n ~ 2, . . . , 9, but this strongly suggests that it is correct. 



22 



Substituting the expression (|102|) into (101), we find that 



Cfc = unless q^^ = . (103) 
This means that most of the Ck vanish. Define an integer a by 



a = N/2 , mod n if is even 
a = {N — n)/2 , mod n ii N,n are odd , 



(104) 



so that in either case < q < n. Then is non-zero only when k 
a,a + n,a + 2n, .... It follows that there are at most 



A/" 



N-a 



n 



(105) 



non-zero vectors in the expansion (|102| ). Here [x\ denotes the integer part 
of X. 

We can write down an explicit, if unwieldy, expression for the elements of 
Cfc by performing a binomial expansion on each product in (|100| ): 



Cfc (3^1 , • • • , 



mi,...,m„ r=l 



(2r-l)mr I ^T+l — Xr —1 



(106) 



where k must satisfy the restriction ( |103| ) and the summation is over all inte- 
gers mi, . . . , nin such that nii + • • • + m.„ = k and 



< rrir < Xr+i — Xr — 1 , for r = 1, . . . , n . 



5.0.4 Numerical tests 

It is not obvious whether these vector are in fact linearly independent. The 
author knows of no reason to suppose they are not, but as a check we have 
numerically calculated the vectors for n = 2 and 3. We can distinguish three 
cases: 

i) n = 2 , A = iTT/2, q = —i, = —1, 

ii) n = 3 , A = i7r/3, q = e-"^""'/^, = 1, 

iii) n = 3 , A = 2m/3, q = e"^^!'^, = -1. 

For all these cases g^" = 1, but only for the second is = 1. Thus 
must be even for cases (i) and (iii), but may be either even or odd for case 
(ii). 

We present the results in Table 1, for = 2, . . . , 16. In every case we 
calculated the vector g with elements (|100| ) for 12 randomly chosen values 
of z\ and then determined (to 15 digit precision) the rank f of the matrix 
with these 12 column vectors. We then numerically verified that each was an 
eigenvector of the six- vertex model transfer matrix T, with eigenvalue (|96[). 
Finally, we calculated the column nullity n of T — AI in the subspace with n 
down arrows. This is the degeneracy of A. In every case we found r = h = Af, 
where TV is given ( |105| ) . Thus at least for these cases V is indeed of dimension 
( |105D and contains all eigenstates. 

We can compare these results with those of Fabricius and McCoy p!5[ |. 
Their 7 is related to our A by A = ^7, so A = — cos 7 and q = —e^^. From 
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N 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


case (i) 







2 




2 




4 




4 




6 




6 




8 


case (ii) 





1 





1 


2 


1 


2 


3 


2 


3 


4 


3 


4 


5 


4 


case(iii) 












2 




2 




2 




4 




4 




4 



Table 1: The dimensions of the space V for n ~ 2 and n = 3, as calculated numerically (a value 
of zero implies that their are no eigenvectors and that (p3|) is not an eigenvalue) . In every case the 



result agrees with ( |105| ) and with the calculated nullity of T — A/, implying that the vectors (106) 
are indeed linearly independent, and that V is the complete eigenspace. 



Table 2 of their paper, when 7 = vr/2, A = 0, n = 2 and = 16, there 
are eight single strings in the = N/2 — n = 6 sub-space, all corresponding 
to the same eigenvalue. This agrees with the above derivation : a = , so 
Af = 8 and the summand in ( [102[ ) is non-zero only when k takes the eight 
values 0,2,4,6,8,10,12,14. 

Also, from Table 7 of the same paper, when 7 = vr/S, A = —1/2, n = 3 and 
N = 16, there are four equal-eigenvalue single strings in the 5^ = N/2 — n = 5 
sub-space. This also agrees with the above, and with case (ii) in the table: 
a = 2, A/" = 4, and /c = 2, 5, 8, 11. 

Some of these results have also been obtained by Braak and Andrei |35], 
who refer to the freedom in the choice of the string centres as "transparent 
excitations". Their table 1 is line 2 of our Table 1 above. 



5.0.5 Conclusions 

Let V be the space that spanned by the non-zero vectors Ck- Then as zi varies, 
the eigenvector g{zi) traces a curve within this space. Each vector Cfc can be 
written as a sum of eigenvectors g{zi), so is itself an eigenvector. V is an 
eigenspace corresponding to the eigenvalue (p^). 

It appears that the vectors Ck are linearly independent and span the full 
eigenspace (within the sub-space = N/2 — n of n down arrows). If so, then 
V is of dimension M. The eigenvalue has degeneracy M, and we have used 
the Bethe ansatz to construct the full eigenspace. 



5.0.6 The case when the string parameters vi,. . . ,Vm are infi- 
nite 

If q satisfies both the restrictions (^) and 

q'' = l , (107) 

then k = and k = N — n in ( |102| ) both correspond to non-zero vectors c^. 
These Ck are the values of g when zi = and zi = 00, i.e. when vi, . . . ,Vm = 
—00 and vi, . . . ,Vm = +00. We have a string of infinite vi, . . . , Vm- 
From ( |100D , we readily find that 

Co(xi,...,X^) = g-i+-2 + -+x™ ^ 
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CM-m{xi.....Xm) = q^-^ q-^'-^—r^ (108) 
These are particular cases of the special vectors reported in equations (16) - 



(22) of in fS^l, and in the equation (|55| ) above. 



Multiple complete strings: the function Q{v) 



We also calculated the null space of the matrix B in (42), algebraically using 
Mathematica. We fixed the values of A and q as in cases (i), (ii), (iii). We 
then allowed n (the number of down arrows) in ( |96|) and (|29|) - (p^ ) to take 
all values from to N ^ for = 2, 3, . . . , 9. Since A.{v) is given by (p6D, we 
could immediately calculate to, . . . ,tN' and form B. We did indeed find that 
the column nullity of B was sometimes greater than one. For cases (ii) and 
(iii), where — 3iA is an integer, we found that the nullity was zero unless n 
was a multiple of 3, meaning that A{v) was not an eigenvalue . When n was 
a multiple of 3 the nullity was (n + 3)/3 and (^) was satisfied provided only 
that Qj was zero when j is not a multiply of 3. Thus 

e"^/^(5(f ) = arbitrary polynomial in of degree n/3 , 

where z = . Whatever the choices of the coefficients of this polynomial, 
it can be factored into n/3 polynomials of degree 1 in z^, each of which is a 
complete string. 

So Q{v) factors into a product of n/3 complete strings of length 3. The 
" centre" (the average value of the three vjs in the string) of each string is 
undetermined. Any such function Q{v) satisfies (|4l|). 

We found corresponding behaviour for case (i), when A = i7r/2: for n even, 
Q{v) factors into a product of n/2 undetermined complete strings of length 
2. There no solutions for n odd. 

We have not attempted to generalize the derivation of this section of the 
eigenvector g to such multiple strings, but presume that it can be done, and 
that one would find the remarkable binomial pattern of degeneracies reported 
by Fabricius and McCoy in their tables 2 and 7 for 'maximum = 8'.|15] 



In the 2n = N sub-space one can impose the "sum-rule" constraint ( 151 ) 
on f 1 + • • • + Vn- (For non-degenerate eigenvalues, with no strings, this will 
automatically be satisfied. If there are strings, it is not necessary, but may 
be convenient, and will still give a complete set of eigenvectors.) If f i, . . . , Vn 
contain only one complete string, then this condition can be used to fix its 
centre. For n = 2 or 3, and = 2n, this gives z" = (— 1)*" exp(— n^A). This 
r is for states symmetric under arrow reversal, 1 for antisymmetric states. 
We have verified numerically that the resulting eigenvectors (|100| ) do indeed 
have this (anti-)symmetry. 



6 Some of the vi, . . . ,Vn equal 

Suppose that all the Sij are non-zero. Then the A(pi, . . . ,pn) are given by 
(^) . Substitute this into (||) and for the moment regard vi , . . . , f „ as arbitrary 
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parameters and xi,...,Xn as fixed. The result is an entire anti-symmetric 
periodic function of ui , . . . , f „ . It must therefore contain the factor 

n smh[{vi-Vj)/2] . (109) 

l<i<j<n 

In principle one can divide this factor out (it is the same for all xi, . . . ,Xn, 
so is just a normalization factor for the vector g). The result is a symmetric 
function of wi, . . . For instance, if n = 2 and A'^ = 4, as in (^) we can 
write Q{v) as a Laurent polynomial m. z = e": 

Q{v) = z-'^{dz^ +ez + f) . (110) 

We can then follow this procedure so as to write all the g{xi,X2) as multi- 
nomials in d,e, f. Define 

A' = 2A = g + g-i , 

^1 = q'^d -e + qf , ^2 = qd-e + q~'^f , 
6 = (d + f)A' - 2e , 
^4 = -4e(d + /) + {d^ + 6df + f + e^)A' - dfA'^ , 
^5 = 2e(e2 - 3d^ - Wdf - 3f ) + {d+ f)A'{d^ + Udf + f + 3e^)- 
eA'\2df + e^) - 3df{d + /)A'^ + de/A'^ , 
then we find that we can normalize the vector g so that 

5(1,2) = e?6 , 5(1,3) = 6^4 , 

5(1,4) =^5 , 5(2,3) =^^2^3 , (111) 

5(2,4) =6^4 , 5(3,4) =e2'e3 . 

Each element is a multinomial of degree 3 in the coefficients d, e, /. Such 
expressions remove the difficulty that occurs when two or more of the vj be- 
come equal: in that case the anti-symmetric factor mentioned above vanishes, 
as does the right-hand side of (|5|), but the above expressions do not. 

This procedure is equivalent to taking the limit of the ratios of the eigen- 
vector elements 5(2:1, . . . ,Xn) in (^) as the vj approach one another. It also 
has the conceptual advantage of making it clear that one does not necessar- 
ily have to calculate the zeros vi, . . . ,Vn of Q{v). Instead one can imagine 
solving the bilinear equations given in (|2^ ) - ([42| ) for the coefficients of the 
Laurent polynomial expansions in of A(v) and Q{v), then using the results 
in equations such as the above. 

Of course, in practice we do not have useful explicit results for the general- 
izations of ( |111| ) to arbitrary n and N, so for anything other than small n, N 
one is forced to use (|5|) and (|43| ) directly. However, from this point of view 
there is nothing remarkable about two of the Vj coinciding: it merely means 
that Q{v) has a repeated zero. If all one wants is the eigenvalue A(f), then 
the alternate form form ( [29| ) - ( p2[ ) of Bethe's equations can be used directly 
as written. 
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7 Extension to the eight-vertex model 



The zero-field eight-vertex model was solved in 1971 by the author[28, by 
extending the functional relation ( ^T|) (with H = 0) to the eight-vertex model, 
provided that 

N = even . (112) 
This restriction applies throughout this section. 

This functional relation method gives the eigenvalues A(f), but not the 
eigenvectors. In 1972, while at Stony Brook, the author derived equations for 



the eigenvectors of the eight- vertex model in a sequence of three papers [M, 37, 



35]. The basic technique was to convert the eight-vertex model to an ice-type 
solid-on-solid model, and then to solve this by an appropriately generalized 
Bethe ansatz. 



Here we use the notation of |31, 37, 38 1 and prefix the equations by I, II 



III, according to in which of the three papers it appears. Papers II and III 
consider the case when the parameter ry satisfies the "root of unity" condition 
(1.9), (II.6.8) or (III.1.9), i.e. 

Lr] = 2miK + im2K' , (113) 

where L,mi,m2 are integers. This is analogous to (^). This restriction is 
not needed in section 6 of I, because the condition (8) therein is sufficient to 
ensure the required cyclic boundary condition from column N to column 1. 

Papers II and III further consider the case when there are integers n, z^' 
such that 

N -2fi = Lv' . (114) 
Modified elliptic theta functions are introduced: 

H{u) = Hjb{u) exp[iTrm2{u- Kf/{4KLrj)] , 

Q{u) = Qjbiu) exp[i7rm2{u- Kf/{4KLri)] , (115) 



H jb{u) , Q jb{u) being the usual Jacobi theta functions (eqn. 15.1.5 of 
These modified functions are periodic of period 2Lrj. The zero-field eight- 
vertex model Boltzmann weights are then given by (1.8) and (II.6.1): 

a = pQ{-2r])@{r] -v)H{r] + v) , 

b = -pe(-2r])H{r]-v)@{r] + v) , 

c = -pH{-27])Q{r]-v)eirj + v) , (116) 

d = pH{-2r])H{r] -v)H{r] + v) . 

One also uses the functions 

h{u) = -h{-u) = H{u)Q{-u) , (117) 



which satisfy 



cl>{u) = [pQ{Q)h{u)f , (118) 
/i(n + Lr/) = (-l)'^i("^2+i)/^(^) ^ (/.(-n) = (/<(n) , (119) 
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h{u + iK') = _e-i'^™i(2«+ii^')/i'7/j(^) ^ 

h{u + 2K) = _e2i'^'»2(«+^)A^/i(u) . (120) 



This work was re-derived by Takhtadzhan and Faddeev (|^9[) using the 
"Quantum Inverse Scattering Method" (QISM). Then in 1982 the author pre- 
sented the functional relation method in sections (10.5) and (10.6) of his 



book[27|. An explicit construction (for even) of the matrix Q{v) for the 
eight-vertex model is given in section 6 of |^l[, and in section (10.5) of p7| . 

The notation in [^] is slightly different from that in I, II, III and ( [|39| ) . 
The restriction ( |113| ) is not made, H(u), Q{u) are the standard theta functions, 
the elliptic integrals K,K' are written as /, and if we write X,v therein as 
\b,vb, thenQ 

XB = 2i{K-r]) , VB = 2i{v-K) . (121) 

As always, there are errors and inconsistencies. This is a good opportunity 
to correct two of them. 

One is a simple but significant typographical error: j + 1 in eqn. (10.5.8) 
should be J — 1, so that it should read 

Sj = s + X{ai + ■■■ + aj^i) . (122) 

Equation (10.5.21) is then consistent with (I.78).[3 

The other is an omission (or at least an over-simplification) by the author 
in the Bethe ansatz derivation of eigenvectors in papers I - III. In (7) of |28| 



each eigenvalue Q{v) of the matrix Q{v) is taken to be simply a product of 
elliptic theta functions. This is corrected in equation (6.10) of p^, and in 



(10.6.8) of where an exponential factor is also included, so that 



N/2 

QBiv) = e^'^^l[Hjbiv-u,)ejbiv-Uj) , (123) 
i=i 

where the suffix B is inserted to distinguish this function from that of III, r 
and ui, . . . ,U]\f/2 satisfy 

T = tt{s - 1 + N + ip) /8K , (124) 

UI + --- + un/2 = {fs-l + 4p)K/2 -i{s-l + N + A^)K' . (125) 

Here are integers, and f = ±1 is the eigenvalue of the operator R that 
reverses all arrows (or spins), and s = ±1 depending on whether the number 
of down arrows is even or odd (and S is the diagonal matrix with entries s). 



We have used the notation of (10.6.7) - (10.6.8) of [27|, except that we have 
converted from the v = vb therein to the present v,Uj (which are those of 
papers I - III) by (|l2l|) . 



^^This still leaves d with a different sign in |27|] from that in I, II, III, but since vertices 7 and 
8 are sinks and sources of arrows, changing the sign of d does not affect the partition function or 
the eigenvalues of the transfer matrix. 

Apart from ±i factors that presumably arise because of the negation of d. 
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For h = N/2, the eigenvalue equation (III. 1.21) should be the same as 
(10.6.1) of After converting to the notation of III, we find that it is the 
same, and is consistent with the other equations (III. 1.1) - (III. 1.23), if extra 
CO factors are included in (III. 1.14), (III. 1.21), (III. 1.23) to make them become 

* = E E ^' /(^i^i' • • • V'C^i, • • • > , (126) 

1=1 X 

Hv-Uj) h{v-Uj) 



, -2 JNki 
U) 



^iNk, = _ Jl h{uj -Um + 2rj)/h{uj -Um- 2r]) . (128) 

m=l 

Here 

uj = e^™/-^ , (129) 
where the integer rh is given by 

2m = mi{s - 1 + N + Ap) +m2{rs - 1 + 4p) . (130) 

The other equations amongst (III. 1.1) - (III. 1.22) remain unaffected, in 
particular e^^^ is defined by (III. 1.17): 

e'^^ = h{uj + v)/h{uj - V) ■ (131) 
We can still define a function Q{v) by (III. 1.24): 

h 

Qiv) = l[hiv-Uj) , (132) 



but it differs from Qb{v) via the exponential factors in ( |115| ) and ( |123| ). Then, 
exhibiting the dependence of the transfer matrix eigenvalue A on f , ( |127| ) can 
be written as 

A{v)Q{v) = LO(p{v-r])Q{v + 27]) +u;-^(j){v + ri)Q{v-27j) , (133) 

which equation replaces (III. 1.25). 

These are in fact the correct equations, not just for h = N/2, but for all 
h satisfying (|114| ).F^ We refer to equations (III. 1.1) - (III. 1.25), with the re- 
placements ( 126D - (|133| ), as the corrected eight vertex Bethe ansatz equations. 

There was an omission in the derivation in paper III. Equation (III. 2. 2) is 
correct as written, but in (III. 3.1) the author should have allowed the more 
general ansatz of including a factor in the rhs, where lo^ = 1. This is 
equivalent to associating this factor with f{l\xi, . . . Xn) and multiplying the 
first term on the rhs of (III. 2. 2) by uj, the second by io~^, and (via the "wanted 



-'^'^ Except that the restrictions ( |125| ), ( |130D apply only for h ~ N/2, which can be regarded as 
the generic case, and not necessarily even then if complete strings are present: we return to this 
point later in the section. 
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terms" and the "unwanted boundary terms" of section 3 of III) to the intro- 
duction of the UJ factors in ( |126 ) - ( |133| ). 



Precisely these a; factors are included in the n = equations (II. 5. 6) and 
(II. 5. 7), where gi = cj}{v + t]) and §2 = 4>{v — rj). They are also included in the 



work of Takhtadzhan and Faddeev |39] 
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As Takhtadzhan and Faddeev comment, [39, after (5.26)] the original equa- 
tions of paper III appear to apply only for the case when m = 0. However, 



it is better than that. One can verify, for all values of n satisfying ( 114 ), 
using ( |12[1| ), that the corrected eight vertex Bethe ansatz equations and (130) 
are unaffected (apart from normalization factors that merely renormalize the 
eigenvector ^), by the following simultaneous substitutions: 

fh ^ fh — 2mi , p' ^ p' — 1 . 

Similarly, if ui is incremented by 2K, then m, p are incremented by 2m2, 
1, respectively. 

The same remarks apply if ui is replaced by any of the ui, . . . , Un- 
We can use this freedom, which is simply the choice of period parallelo- 
grams for the zeros ui, . . . ^Un, to increment fh by any even integer. If L is 
odd, this means that we can construct any choice of u: by such shifts; while 
if L is even, we can construct half the choices. So the 1973 papers do cover 
approximately three-quarters of the cases! 

More recently, related equations for the eigenvectors of the eight-vertex 



model have been studied by Felder and Varchenko|41], and by Deguchi 



7.0.7 Alternative form of Bethe's equations 

From ( 120[ ) and ( |132| ), remembering that is even, 



Q{v + u + iK') = (-l)"e-™i(2''+*^')/^''Q(?; + u) , 

Q{V + U + 2K) = (-l)"e2™2{«.+i^)/Lr,g(^^^) ^ 

(l){v + iK') = e-*^^™i(2^+*^')/i'7 0(^) ^ (134) 

where 

u = {Ul^ h Un)/h . (135) 



From (116), or from ( |114D and ( |133| ), the function A{v) satisfies the same 
quasi-periodicity relations as those above for (j){v). 

The functions Q{v), (l){v),A{v) are all entire. Using the v ^ v + 2K quasi- 
periodicities above, it follows that there exist coefficients qj,fj,tj such that 

j 
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Our LO is e^^TTzm/Q -^^ 

the notation of Takhtadzhan and Faddeev, and it seems that the I of our 



equation (126) must be the —I of their equation (5.23). 
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j 

where j takes all positive and negative integer values and j' = j if n is even, 
while j' = j— 1/2 if h is odd. 

Now using the v v + iK' quasi-periodicities, we find that the coefficients 
in these series must satisfy 

If gO) • • • Qh-i are known, then the simple periodicity relation ( |137| ) deter- 
mines all the other qj. Similarly, all the tj are determined by to, . . . ,tN-i, 
and the known by /o, . . . , /n-i- 

These series are convergent for all finite v. Substituting them into ( 133 ) 
and equating coefficients, we obtain 

oo 

J2 Bj,mqm = (138) 
m=— oo 

for all integers j, where 

Tj _j. j_ iiTm2(8n—N)r]/2LK\~ r _L f 1 

^j,m — '^j—m ~r c ['-^Jj~m+m2i'' ' ^ J j—m—m2u' \ j 

where 

- in(m+2m'—j)ri/K —2innm2u/LK 

and m' = m if n is even, m' = m— 1/2 if n is odd. 
Define 

= e-^"^'/2ni^^. , (139) 

Bj^^ = ^_iy+m^nj''K'/2(N+n)K ^-nm'^K'/2nK j^^^^ ^ (14q) 

Then ( |138[) becomes 

oo 

{-irB^^mlm = (141) 

m=— oo 

and g^, -Bj,m satisfy the periodicity relations 

qm+h='^m 1 Bj+N+fi,m+h = Bj,m ^ m,j . (142) 

We can therefore define discrete Fourier transforms ija, Bj^a such that 

n— 1 oo 
a=0 m=—oo 

m, a being integers in the ranges — oo < m < oo and < a < n. Then (141) 
becomes 

n-l 

= • (144) 

a=0 
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Here j can take any integer value, but from (142) 



Bj+N+h,a = Bj^a , (145) 

SO there is no loss of information in restricting j in ( |144D to lie in the range 

< j < N + h . 

Hence (|144| ) can be regarded as a set of N+n homogeneous linear equations 
for the n unknowns qa- It is also linear in the coefficients tj, and all of these can 
be determined immediately and linearly from ( |137[) if we know to, . . . ,tAr_i. 
Thus to, ... , t]\f-i play the role of generalized eigenvalues, as to, ■ ■ ■ ,tj\f do for 
the six- vertex model in (^), (0)- Since only the ratios of the Qa enter the 
equation, we have N + h equations for a total of + n — 1 unknowns. Unlike 
the six-vertex model, this set is over-determined, presumably because of the 
quasi-periodicity constraints satisified by the elliptic functions. They must of 
course have solutions. 

Many of the remarks made following equation (^8|) for the six- vertex model 
extend to the eight-vertex model. One can construct an -|- n by n matrix 



B with elements -Bj,a, and write ( 144 ) as Bq = 0. Thus B must have rank at 
most n — 1, and q is the column null vector of B.|^ If the eigenvalue A(u) is 
degenerate (for all v), then the rank of B will be less than h — 1. There will 
then be more than one solution for q, and hence of ( |133| ) for Q{v). One can 
expect such behaviour in the situation that we shall now discuss, i.e. when 
Q{v) contains one or more complete strings. 

7.0.8 Similarities to the six-vertex model: strings 

The eight-vertex model is a generalization of the zero-field six-vertex model, 
and many of the remarks we have made about the six-vertex model continue 
to apply. There are still very special eigenvectors with eigenvalues of the form 
(^), namely the eigenvectors discussed in papers I and II, corresponding to 



n = in (127) . For N even, explicit expressions for the matrix Q are given 
in ||2^ and in section 10.5 of |^^. From these it follows that each eigenvalue 
Q{v) is a product of N/2 elliptic h{u) functions, so the only way we can get 
these simple eigenvalues is for all the ui, . . . to be grouped into complete 
strings. These are of length M, where 

M = L/2 if L is even , 

= L if L is odd . (146) 

Each string consists of M zeros, say ui, . . . ,um, differing sequentially by 2r], 
i.e. for j = 1,...,M 

Uj+i = Uj + 277 , (147) 
interpreting um+\ as u\, modulo 2Mri. 



^^It is a compication that B depends on u. This disappears if m2 — 0: possibly it can be removed 
for other values by using the periodicity in integer multiples of r\ (rather than as the basis for 
the discrete Fourier transforms. 
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In fact any eigenvector and eigenvalue given by (III.l.l) - (III. 1.23)[] with 
h 7^ N/2, must also be present in the case n = N/2, differing from it only in 
the subtraction or addition of complete strings. The eigenvectors for h ^ N/2 
lie in the eigenspace of any other allowed value of n which is closer to (or 
equal to) N/2, at least for n < N/2. 

For any values of L, mi, 1712, the condition ( |114| ) is satisfied by n = N/2, 



so this is the generic case. The eigenvalue equations (|127|) - (131) then apply 



for all r]. The eigenvector equations, notably ( |126|) , depend on i] satisfying 
( |113D , but one can approach arbitrarily close to any desired value by taking 
L, mi, 1712 sufficiently large 

Other values of n only occur if ( |113| ) is satisfied for L not greater than A^. 

Again, there are technical difficulties about handling Bethe's equations 
when there are complete strings. The remarks of section 4 extend from the 
six-vertex to the eight-vertex model. One should replace ( |128D by the two 
equations 



uj e 



2JNk 



n i-tm,j/tj,m) , (148) 
m=l,mj^j 

tj,m h{Uj -Um + 2r]) = tm,j h{Um - Uj + 2r]) , (149) 
and write (III. 1.20) as 

A{P) = ep n tP^^^Pj ' (150) 

l<j <m,<h 

where the renormalization factor C is the same for all permutations P and is 
the product of selected tj^m factors, one from each string, for which hium — 
Uj + 2r]) vanishes. 

If Uj does not belong to any string, then in ( |148[ ) we can take tj^rn = 
h{um — Uj + 2r]). The contribution to the rhs from the Um that do lie within 
strings cancels out, leaving a reduced equation where j, m only take non- 
string values. If Uj does belong to a string, say to ( |147] ) , then ^12,^23; • • • ,tM,i 
vanish but their ratios remain finite. If one fixes one of the Uj within the 



string, then the rest are determined, the e^^^ are given by ( |131| ), and the 
ratios of ti2, ^23, • • • , ^M,! are determined by ( |148|) for j = 1, . . . , M. From 
(150), some of the coefficients A{P) involve ti2,t23, ■ ■ ■ ,tM,i only via these 



ratios, so are finite and non-zero. The other A{P) are of linear or higher 
order in ti2, ^23, ■ ■ ■ , tu,!, so vanish. 

Again, one is free to vary each of the string centres at will. As one varies 
these parameters, and the disposable parameters s, t in III, the eigenvector 
^ will trace out a surface <S in the eigenspace of the eigenvalue A. If the 
eigenvalue is unique, these variations will merely change the normalization. 
If it is degenerate, S will lie in the eigenspace appropriate to this value of n, 
and we expect the vectors on S to span this eigenspace. 



^^including the modifi catio ns ( 126 ) - (|l2§| ) 

^^The sum over I in (|l26| ) is a discrete Fourier transform: for general values of 77 it may 



appropriate to replace it by a continuous one. 
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Since the explicit construction of Q{v) given in section 10.6 of gives 
n = N/2, we expect this to be the generic case, giving all eigenvalues and a 
complete set of eigenvectors of T{v). For h satisfying (|114|) , but not equal to 
N/2, one expects to only observe the degenerate eigenvalues, and to obtain 
only a sub-space of the eigenspace of each. Recent numerical results support 
these expectations. [^] 

Of course, one may have particular reasons for fixing the string centres at 
particular values, as Fabricius and McCoy did for the six- vertex model. 1 14 1 
[17| An obvious choice (for h = N/2) is to fix them so that Q{v) is the 



eigenvalue of the matrix Q{v) constructed in section 6 of [29| and in section 
(10.5) of [p^ j. Equivalently, one can require that they be fixed to their limiting 
values (again, for h = N/2) as rj approaches the "root of unity" value ( 113 ). 
However, these considerations lie outside the Bethe ansatz for a fixed value of 
tj. The Bethe ansatz is complete without them: the arbitrariness in the string 
centres (and in s, t) is a reflection of the degeneracy of the eigenvalues of the 
tranfer matrix, and the resulting non- uniqueness of the eigenvectors. 



7.0.9 The six-vertex model limit 

In the limit when the elliptic modulus k (or the nome q) goes to 0, 1 or 
oo, the elliptic functions become trigonometric functions and the eight-vertex 
model becomes the zero-field six-vertex model. Much of the working of this 
section can be adapted at once to the six-vertex model, except that some 
of the ui, . . . ,uj^/2 may become infinite. In this way the resulting six- vertex 
model function Q{v) can have any number n < N/2 of finite zeros, and the 
oj factors in (|127] ) can be related to those in (|73|). It should be possible to 
obtain all the six-model eigenvalues from the those of the eight-vertex model 
by taking such a limit. 

There is a problem with the eigenvectors. For n ^ N/2, the zero-field six- 
vertex model eigenvalues occur in degenerate pairs, one in the sub-space with 
n down arrows, the other in the arrow-inverted sub-space with N — n down 
arrows. Thus two eight-vertex model eigenvalues, with opposite spin-reversal 
symmetry, must coalesce. Their sum and difference will then be the six-vertex 
model eigenvectors in the two sub-spaces. Only in the n = N/2 sub-space can 
one expect to obtain the six-vertex model eigenvectors directly as limits of 
those of the eight-vertex model. 



7.0.10 The sum rule 



The constraint (125) applies to the eight- vertex Q{v) functions obtained by 
the explicit construction in section (10.6) of (2^], which have h = N/2 zeros. 
If the eigenvalue A.){v) is non-degenerate, i.e. if there are no exact complete 
strings, then Q{v) is uniquely defined by Beth's equations, so h will be equal 
to N/2, and ( |125| ) will automatically be satisfied. 

However, if there are strings, even if n = N/2, then one can shift the string 
centres arbitrarily and Bethe's equations will be unaffected. The resulting 
eigenvector will not necessarily be an eigenvector of R and S, and ( [L25D will 
not in general be satisfied, but it will nevertheless be a valid eigenvector of 
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the transfer matrix T{v). 

So we conclude that, for all rj satisfying ( [L13D , a complete set of 2^ eigen- 
values and eigenvectors can be obtained by taking h = N/2 and observing 
the constraint ( |125| ) (though even then the eigenvectors will not necessarily 
also be eigenvectors of R and S for all values of the string centres and the 
disposable parameters s and t |p3|). Further eigenvectors can be obtained 



by abandoning these constraints, while of course retaining (114), but these 
eigenvectors will lie in the eigenspaces obtained with the constraints, so do 
not extend these eigenspaces. 

These remarks extend to the six- vertex model limit. In the n = N/2 sub- 
space, non-degenerate eigenvalues must satisfy the analogue of ( |125| ), which , 
using (|121| ), is 

VI + --- + VN/2 = iTT{rs -1 + N + 4p)/2 . (151) 

These f i, . . . , V]\f/2 are those of sections 1 through 6, p is an arbitrary integer, 
r the eigenvalue (±1) of the spin reversal operator, and s = (—1)^/^ the 
arrow parity of this state, with N/2 down arrows. This is equation (17) 



of ||4j], and equation (1.44) of ||17|]. Again, if the eigenvalue is degenerate 
and ^^i, . . . ,f7v/2 contain one or mpore complete strings, then (|151| ) will not 
necessarily be satisfied. I am indebted to Barry McCoy and Klaus Fabricius 
for correspondence on this and related matters concerning the Bethe ansatz. 



7.0.11 Other possible difficulties: bound pairs 

Because of the double periodicities of elliptic functions, the technical prob- 
lems in the six- vertex model associated with zeros vi, . . . ,Vn going to infinity 
cannot occur in the eight- vertex model: the corresponding ui, . . . ,Un can be 
restricted to a period parallelogram. There appears to be no "beyond the 
equator" problem. All the states are accounted for by taking h = N/2. The 
equations may well have solutions for h > N/2: this would correspond to 
adding complete strings to Q{v). 

One could still have bound pairs analogous to (^), when — ui = U2 = rj, 
^iki ^ g-ifc2 ^ , e^('=i+^2) ^ _i ^ (152) 

and the ratio t2i/ii2 vanishes. One would expect to handle this in the 
same way as for the six- vertex model, using ( |148| ) for j = 1,2 to calculate 
e'^^^^ti2/t2i and e'^^^^t2i/ti2, then substituting these into the suitably renor- 
malized equations for A(P) and the eigenvector ^. 

Although we have not observed the phenomenon, it is conceivable that 
two or more of ui,. . . ,Un, say ui, . . . ,Up, could coincide at some arbitrary 
value. This would require either generalizing the argument of section 6, or 
first dividing (III. 1.16) by the product of h{ui — Uj) over 1 < i < j < p, and 
then taking the limit where ui, . . . ,Up become equal. 



8 Summary 

We have presented the coordinate Bethe ansatz equations with some care, 
trying to avoid (or at least signpost) the problems that occur when some of 
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the variables are zero or infinity. Perhaps the essential point of this paper is 
that for the six- vertex model the Bethe ansatz equations are (|5|) - (p!7[). It 
seems that one can always choose the coefficients A{P) to be given by (^3|), 
but this is not necessary if enough of the Sij vanish (as they do when the Vj 
are equal and infinite). All that is necessary for g to be an eigenvector is that 
® " (0) {^^ their appropriately renormalized forms) be satisfied, with g ^ 0. 

For the six-vertex model, in section 3 we have discussed the situation that 
arises when some of the Bethe zeros Vj are infinite, and how this leads to a 
reduced Bethe equation containing w factors. In section 4 we show that this 
is the key to resolving the "beyond the equator" problem and to constructing 
the Bethe eigenvector for n > N/2. We also show how to cope with the 
problem of a bound pair, when two of the momenta are infinite but their sum 
is finite (equal to an odd integer multiple of in). 

In section 5 we look at the problem discussed by Fabricius and McCoy, 
when there are one or more exact complete strings. We show that the vj are no 
longer uniquely determined, because the string centres can be chosen at will. 
Nevertheless, the Bethe ansatz equations are satisfied, and we show how to 
construct the (necessarily non-zero) eigenvector by working with appropriate 
ratios of the vanishing tij . 

In particular, we have found the solutions of the Bethe ansatz correspond- 
ing to all the vi,. . . ,Vn lying on a single complete string. The Bethe ansatz 
equations do not define the string centre (the average of wi, . . . f„). This is to 
be expected: it is a direct consequence of the eigenvalue A being degenerate, 
which means that there is more than one eigenvector, and hence more than 
one solution of the Bethe ansatz. We show that the ansatz can be used to 
construct a complete set of eigenvectors, spanning the eigenspace. 

One should of course go on to study more complicated situations, where 
there are more than one complete strings, and not all the vj belong to a string. 
We expect the above methods to generalize to such cases, but the algebra may 
well be complicated. 

In all the cases dealt with in sections 3 to 5, one can always use (^) (with an 
appropriate choice of the normalization factor C in (^)) as written, each term 
in the summand being finite and the sum being finite and non-zero. The only 
"limit" is that of recognizing that one has the set of rational equations (20), 
(^) - ( ^6[ ) to solve for the variables e^^ , e'^^^ and tij . Some of these variables 
may be zero or infinite, while what one wants is their ratios or some other 
product of powers, which are finite and non-zero. These ratios or products can 
of course themselves be regarded as variables. This is only a generalization 
of the usual practice of including the "point at infinity" in the domain of the 
variables. 

In section 6 we touch on another problem, namely what happens if two 
or more of the Vj are equal and finite. (The question of what happens when 
they are equal and infinite is different, actually easier to resolve, and is dealt 
with in section 3.) In general this case really does seem to demand that one 
take a limit in the expression (^) for the eigenvector g, since the terms in the 
summand are finite and non-zero but cancel one another in pairs, so that their 
sum is zero. This is different from and less satisfactory than the other cases, 
but we remark that in fact we never encountered this case in our numerical 
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experiments, and it is not the case discussed in Q - |jT^. It is not clear that 
it ever actuahy occurs. 

In section 7 we discuss the zero-field eight-vertex model with an even 
number of columns, and give the needed corrections to the coordinate Bethe 
ansatz equations of section 1 of |38]. We indicate how the string and infinite 
momenta problems can occur also for this model, and how to resolve them. 

If one wants to make a specific choice of the string centres, particularly if 
one wants to ensure continuity as A or r] passes through the "root of unity" 
value ( pOD or ( |113D , or (equivalently) if one wants Q{v) to be the eigenvalue 
of the matrix Q{v) of (86) of |Q, or (10.5.31) of then clearly one should 
use the results of Fabricius and McCoy. They have addressed this problem in 
a series of well-presented papers, and have systematically exhibited the con- 
nections to the sl2 loop algebras. However, they do use the provocative title 
"Bethe's equation is incomplete .-.".[ISl If all one wants to do is to diagonalize 
the transfer matrix (or the XXZ hamiltonian) , obtaining all the eigenvalues, 
their degeneracies and eigenspaces, then it seems that there is no need to look 
further than the Bethe ansatz. At least for the cases studied in this paper, 
the Bethe ansatz is in fact complete. 

The alternative forms (41), ( |138D of Bethe's equations are themselves gen- 
eralized eigenvalue equations in which the n + 1 or n independent coefficents 
Qj are the elements of the eigenvector. They have some advantages over (44), 
(H), (HH), ( |149D , being hnear in the qj and the tj. The "infinite vj" prob- 
lem merely corresponds to some of the coefficients qj vanishing.^ It seems 
surprising, but while the author did not believe these equations to be new, 
he has been unable to find any previous paper where they have been written 
down. 



Appendix A 

Suppose for the moment that ^(1, . . . ,n) is non-zero. Let Omj be A{Pfn,j), 
where Pm,j is the permutation where j is removed from its place in the se- 
quence 1,2,... ,n and replaced immediately after m, or immediately before 
m + 1. Thus 

am,j = A{l,...,m,j,m + l,...,j + l,...,n) if m< j , 

am,j = A{l,...,j + l,...,m,j,m + l,...,n) if m > j . (Al) 

In particular, 

ao,j = MjA, • • • ,i - 1,J + 1, • • • ,n) , 

= ajj = A{l,...,n) , 

an,j = A{1,.. . , j - 1, j + 1, . . . ,n,j) . 

Fix j at some value between 1 and n. Then there are n — 1 equations of 
the set that involve only the n distinct coefficients aoj, . . . , a-nj, namely 

But in general one still needs to calculate the zeros vj or Uj in order to obtain the eigenvector. 
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for m = 1, . . . ,n, m ^ j. Also, from (^) we find that 

e'^^'^anj = ao,j . (A3) 

. Together, these give us n linear homogeneous equations in n unknowns. 
Since at least one of the unknowns, namely ctj-ij = ctjj = A{1, . . . ,n), is 
non-zero, the determinant of the matrix of coefficients of these n equations 
must vanish. This determinant is easily obtained, giving 

n n 

e^"^'^ n ^J^rn = (-1)"-' n (A4) 

m=l,m^j m=l,m,jtj 

and this must hold for all the possible values 1, . . . , n of j. 



Permuting the indices l,...,n merely rearranges the n equations (A4), 
so our initial assumption that A(l, . . . ,n) is non-zero is irrelevant: to derive 
( |A4| ) it is sufficient that any one of the coefficients A{pi, . . . be non-zero. 
This must be so for y to be a non-zero vector. 
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